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Abstract 

We revisit a single-server retrial queue with two independent Poisson streams (corresponding to 
two types of customers) and two orbits. The size of each orbit is inhnite. The exponential server 
(with a rate independent of the type of customers) can hold at most one customer at a time and there 
is no waiting room. Upon arrival, if a type i customer (i = 1, 2) hnds a busy server, it will join the 
type i orbit. After an exponential time with a constant (retrial) rate an type i customer attempts 
to get service. This model has been recently studied by Avrachenkov, Nain and Yechiali [3] by solving 
a Riemann-Hilbert boundary value problem. One may notice that, this model is not a random walk 
in the quarter plane. Instead, it can be viewed as a random walk in the quarter plane modulated by 
a two-state Markov chain, or a two-dimensional quasi-birth-and-death (QBD) process. The special 
structure of this chain allows us to deal with the fundamental form corresponding to one state of 
the chain at a time, and therefore it can be studied through a boundary value problem. Inspired by 
this fact, in this paper, we focus on the tail asymptotic behaviour of the stationary joint probability 
distribution of the two orbits with either an idle or busy server by using the kernel method, a different 
one that does not require a full determination of the unknown generating function. To take advantage 
of existing literature results on the kernel method, we identify a censored random walk, which is an 
usual walk in the quarter plane. This technique can also be used for other random walks modulated 
by a finite-state Markov chain with a similar structure property. 

Keywords: Retrial queue • Random walks in the quarter plane • Random walks in the quarter 
plane modulated by a finite-state Markov chain • Censored Markov chain • Stationary distribution • 
Generating function • Kernel method • Exact tail asymptotics 
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1 Introduction 

In this paper, we revisit a single server retrial queue system with two orbits and no waiting room, 
which has been studied by Avrachenkov, Nain and Yechiali [3]. The analysis in [3] is based on the 
solution of a Riemann-Hilbert boundary value problem, while our focus is on exact tail asymptotics for 
the joint stationary distribution of the two orbits under a busy or idle state of the server, using the 
kernel method, a different method that does not require a full determination of the unknown generating 
function. In this system, there are two independent exogenous Poisson streams (representing two types 
of customers) flowing into the server, and the server can hold at most one customer at a time. Upon 
arrival, if a type i customer finds a busy server, it will join its orbit and wait for retrial at a specified 
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exponential rate for the customers of type i. Such a queueing system could serve as a model for two 
competing job streams in a carrier sensing multiple access system, and it has an application in a local 
area computer network (LAN) as explained in [3]. 

Retrial queueing systems have been attracting researchers’ attention for many years (e.g., [Il[2l[5l[25] 
and references therein). Much of the previous work lays the emphasis on performance measures, such as 
the mean size of the orbit, the average number of the customers in the system, the average waiting time 
among others. We also notice that stationary tail asymptotic analysis has recently become one of the 
central research topics for retrial queues due to not only its own importance, but also its applications 
in approximation and performance bounds. For example, in [24], Shang, Liu and Li proved that the 
stationary queue length of the M/G/1 retrial queue has a sub exponential tail if the queue length of 
the corresponding M/G/1 queue has a tail of the same type. Kim, Kim and Kim extended the study 
on the MjGjl retrial queue in [llj by Kim, Kim and Ko to a MAP/G/l retrial queue, and obtained 
tail asymptotics for the queue size distribution in [T2|. By adopting matrix-analytic theory and the 
censoring technique in m. Liu, Wang and Zhao studied the M/M/c retrial queues with non-persistent 
customers and obtained tail asymptotics for the joint stationary distribution of the number of retrial 
customers in the orbit and the number of busy servers. 

Most of the studies on retrial queues assumed a single type of customers flowing into the system, 
and references on retrial systems with multi-class customers are quite limited. The model studied by 
Avrachenkov, Nain and Yechiali in |3] and again in this paper is such a system. This model is an 
example of the two-dimensional QBD process (for example, see Ozawa [22j). or the random walk in the 
quarter plane modulated by a two-state Markov chain (another example of retrial queues having this 
structure is Li and Zhao m)- In [3|, the authors showed how this modulated model is converted to a 
scalar fundamental form, which can be solved in terms of a Riemann-Hilbert boundary value problem 
(BVP) due to its special structure of this system. Motivated by this, we extend their research on this 
model by considering the tail asymptotic behaviour of the stationary joint probability distribution of 
the two orbits with either an idle or a busy server, by using the kernel method — a different method 
that does not require a full determination of the unknown generating function. For more details about 
the kernel method, readers may refer to liiiaiisKiTKiBi. We point out that tail asymptotic properties 
for Markov modulated or more general block-structured random walks have also been studied by using 
other methods, for example in |lll[23l[l3l|20|. 

The main contribution in this paper includes; (1) the characterization of the tail asymptotic prop¬ 
erties in the joint distribution for a large queue i {i = 1,2) with either an idle or a busy server. A total 
of three types of properties are identified (see Theorems 16.2116.31 and 16.41 for the case of a busy server, 
and Theorems 16.51 and 16.61 for the case of a idle server); and (2) an illustration on how to convert a 
matrix-form fundamental form for the Markov modulated random walk into a (scalar) functional form 
corresponding to one state of the chain, through a censored Markov chain or solving the matrix-form 
fundamental form (see remarks in the last section). Therefore, it can be studied by the kernel method. 

The rest of the paper is organized as follows: Section 2 provides the model description; Section 3 
identifies the censored random walk in the quarter plane; dominant singularities of the unknown gener¬ 
ating function are located in Section 4, while the detailed asymptotic property of the unknown function 
at its dominant singularity is discussed in Section 5; exact tail asymptotic properties for stationary prob¬ 
abilities of the system, which are our main results, and their detail proofs are presented in Section 6. 
Concluding remarks are made in the final section. 
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2 Model description 


In this paper, we consider a single server queueing system with two independent Poisson streams of 
arrivals and two retrial orbits, the same system studied in [3]. Following [5], the two arrival rates are 
denoted by Aj, i = 1, 2, with A = Ai + A 2 . The server can hold at most one customer at a time (without 
a waiting room). It means that when the server is busy, an arriving type i customer will join in orbit 
i of infinity capacity. Retrials from all customers in orbit i for service are characterized by a Poisson 
process with constant rate /r*. The service time for each customer is independent of its type and follows 
an exponential distribution with rate //. The retrial mechanism imposed can be a model when only the 
customer at head of the line (orbit) is allowed for retrial. 

Let I{t) be the state of the server (either idle or busy), or the number of customers in the server, 
and let Qi{t) denote the number of customers in orbit i at time t for i = 1,2. Then, it generates a 
continuous time Markov chain X{t) = {{Qi{t),Q 2 {t), I{t)) : t G [0,oo]} on the state space {0,1,.. .} x 
{0,1,. ..} X {0,1}. From Avrachenkov, Nain and Yechiali [3], we know that the system is stable if and 
only if A(Ai + /ri) < and A(A 2 + (^ 2 ) < Under the stability condition, the unique stationary 
probability vector for the system is denoted by Il^.n = , T^m,ni^)) for m, n = 0,1,.... For the 

purpose of hnding the stationary distribution, we consider the corresponding discrete time Markov chain 
through the uniformization technique. Without loss of generality, we assume that A + /r + /ri+/r 2 = l- 
For this discrete time chain, a transition diagram, partitioned according to the state of the server, is 
depicted in Figure 1, where 


. ^ _ 4(1) _ 4(2) _ .(0) _ 
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n — 


m,o 


A 


Ail + Ai2 


3 



Figure 1 Matrix transition diagram 

We define the probability generating function (PGF) P^^\x,y) for the stationary probabilities 
'^m,n{k) as 


|x|<l,|y|<l, A; = 0,1, 


rri,=0 n=0 


and denote 


P{x,y) = {P^^\x,y),P^^\x,y)). 

Following the idea in Fayolle, lasnogorodski and Malyshev [6], we can obtain the (matrix-form) funda¬ 
mental form for the Markov modulated random walk in the quarter plane: 

P{x,y)H{x,y) = P{x,0)Hi{x,y) + P{0,y)H2{x,y) + Uo^oHo{x,y), (2.1) 


where 


Hix,y) = -h{x,y), 

Hi{x, y) = -h{x, y) -F hi{x, y)y, 

H 2 {x, y) = -h{x, y) -F h 2 {x, y)x, 

Ho{x,y) = hQ{x,y)xy+ h{x,y) - hi{x,y)y - h 2 {x,y)x 


with 


1 1 


h{x,y) = xy 




i=-lj=-l 
1 1 


hi(a:, 2 /) = ^ ^ Am ^*2/^ - I) , 

i=-lj=0 
1 1 

h(x,y) = y(YY 

i=0 j=-l 
1 1 

ko(x,y) = YY^U^'^^ 


i=0 j=0 


For detailed derivation, see the work in [9]. 
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Remark 2.1 It is worthwhile to mention that the fundamental form (1.3.6) in [6] is for the generating 
function excluding boundary probabilities, while ours is for the complete joint probability vector. In [9], 
it is pointed out that these two forms (for Markov modulated random walks) are equivalent. In fact, 
for /c = 0,1 let 

OO OO 

m=l n=l 

and TT{x,y) = {x, y), {x, y)), then 

-7r{x, y)h{x,y) = 7r(x, 0)/ii(x, y) + 7r{0,y)h2{x,y) + Uofiho{x,y) 


by noticing that 

OO 

Trf'\x) = 7r^^\x,0) = ^ 7rm,o{k)x 

m=l 

OO 

A^\y) = ^^^^(0,2/) =^7ro,n(A;)2/" 

n=l 


For the retrial queueing system studied in this paper, after some calculations, we have 


H{x,y) 


+ Ml + y2)xy 


=(oo)' 

Hence, the fundamental form (12.ip can be simplified as 

P{x,y)H{x,y) = P{x,0)Hi{x,y) + P{0,y)H2ix,y). 


Equivalently, 


^ ^ ^ \ -y^xy -[X2xy^ + Xix^y - {X +fi)xy] 

= (pW(x,0),pW(x,0)) +{P^^\o,y),P^^Ho,y)) 


or. 


(A + Mi + fi2)P^°Hx,y) - nP^^'>{x,y) = fi2P^^\x,0) + fiiP^°\o,y), (2.2) 

{Xxy + my + H 2 x)P'^°\x,y) + [Aix + A 2 y - {X + n)]xyP’'^\x,y) = H2xP^°Kx,0) + (2-3) 

(12.21) and (j2.3p are identical to equations (18) and (19) in [3], derived from direct calculations. 
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3 Censored Markov chain 


One may notice that equations (1^ and (|2.3p provide a relationship between generating functions 
for an idle server and for a busy server. Therefore, we start our analysis for a busy server since in 
this case, the censored Markov chain can be expressed explicitly. This censored Markov chain is a 
random walk in the quarter plane, which has been extensively studied in the literature. To this end, 
we first consider the uniformized discrete time Markov chain of the continuous time chain X{t) for the 
retrial model with uniformization parameter 6 = A + /r + //i + /Z 2 = 1- We partition the transition 
matrix P of the uniformized chain according to the server state and then consider the censored chain 
to the set of states of a busy server. Specifically, let X{n) = {{Qn\ Qn\ In)} be the uniformized 
chain on the state space {0,1,...} x {0,1,...} x {0,1} and let S = {0,1,...} x {0,1,...} x {1} and 
^^ = {0,l,...}x{0,l,...}x {0}. Partition the transition matrix P according to E and its complement 
into: 

E^ E 

E-fPoo Poi\ 

E \Pio PiiJ^ 

where using the lexicographical order for states of {Qn\Qn^), Pij can be expressed as 


with 


-Poo = 


f A, \ 

Ai 

Ai 

\ '■■) 


Poi = 


f Co 


Pio = 


Co 


V 


/ /^ + + M2 


Ao = 


M + Ml 


M + Mi 


V 


Bo = 


( A \ 

M2 A 
M2 A 


V 


= 


Ml + M2 A 2 


/ Ml 

V 


/^o \ 

Bi Bo 
Bi Bo 

V 


Pll = 


(Do Di 

Do Di 


\ 


= 


/m + M2 \ 
M 

M 


Ml 


V 


Co = 




( M 


M 


Do = 


Ml + M2 A 2 


Ai 

Di=\ Ai 


Notice that Pqo is diagonal, it is straightforward to have the fundamental matrix of Pqo as follows: 


Poo = ^ Poo = diag(^O)^i)^i) ■ ■ ■), 

n=0 
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where 


io = diag 
Ai = diag 


_^ ^ 

A’A + /i2’A + a/2’"7 ’ 


1 


1 


1 


X fll A + /ii + /i2 A + /Ti + /i2 


Furthermore, notice that Pio is also diagonal, and therefore the censored chain to E can be easily 
computed as 


p{E) 


^ Dq + hAqBq Di ^ 

IxAiBi Dq + fiAiBo Di 

^AiBi Dq + ixAiBq Di 

V 


This censored chain is an example, referred to the simple random walk, of the random walks in the 
quarter plane studied in [6], whose transition diagram is depicted in Figure 2. 



Figure 2 Transition diagram of the censored random walk 


In our case. 


„ _ „(i) _ „(2) _ „(0) _ X „ _ „(i) _ J2) _ „(0) _ X . 
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Pofl = + At2 + 


where jli = nfii for i = 1, 2. 

Let a = A + /ri+/i2 = l — h A* = aXi for i = 1,2. It follows from [3] that under the system 
stability condition (for the retrial queue), at least one of Ai < fii and A 2 < /I 2 holds. Without loss 
of generality, we assume that Ai < jli throughout the paper. For this censored random walk, the 
fundamental form (equation (1.3.6) in [6]) is given by: 


- h{x,y)7r^^\x,y) = hi{x,y)7r\^\x) + h2{x,y)7:^^\y) + ho(x, y)7ro,o(l), (3.1) 

where 

1 1 

h{x,y) = xy(^ ^ ^ Pijx^y^ “ ^ “(a;)?/^ + b{x)y + c{x), (3.2) 

i=-ij=-i 

1 1 

hi{x,y) = x(^ ^ '^pf]x'-y^ -1^ = ai(x)y + 6 i(x), (3.3) 

i=-l j=0 
1 1 

h2{x,y) = y(^J2 = a 2 (a;)y^+ 62 ( 3 ;)^+ C 2 (x), (3.4) 

i=0j=-l 

1 1 

ho{x,y) = '^J2ptj^''y^ “ ^ = Mx)y + bo{x), (3.5) 


a{x)=Po,ix, b{x) = p-ifi - {1 - po^o)x + pifix"^, c{x) = po-ix] 
ai{x)=PQ^ix, bi[x) = pI{ q - [1 - pl)J^)x + p\J)X ; 

02(3:) =Po}_, b 2 ix) = - 1 +ri^ii 3 :, C 2 (x) = p^qI^] 

ao{x) =Po% bo(x) = p[°ix - (1 -pJJ’J). 

In [ISKIT], a kernel method has been promoted for studying exact tail asymptotic properties for 
random walks in the quarter plane. In the following, we apply this method to the retrial queue model 
to explicitly (in terms of system parameters) characterize regions on which different tail asymptotic 
properties hold. First, based on the fundamental form in (|3.1I) . asymptotic properties at the dominant 
singularity for the generating function (x) or (y) are obtained, based on which regions of different 
exact tail asymptotic properties for probabilities 7rm,n(l) with a fixed value of n or m are identihed 
through a Tauberian-like theorem (Theorem 16.ip . Then, based on the relationship given in (|2.2p . the 
generating functions (x) and vr® (y) are analyzed, and characterization of the exact tail asymptotic 
properties for 7rm,n(0) is provided. 

4 Dominant singularity of 

Since discussions for dominant singularities of the two functions 7rj^^(x) and ir^^^ly) are repetitive, we 
only provide details for 7r[^^(x). 

According to [13 E], the dominant singularity of 7 rj^^(x) is either a branch point of the Riemann 
surface defined by the kernel equation h(x,y) = 0, or a pole of the function 7rj^^(a:). The following two 
subsections are devoted to these two cases, respectively. 
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4.1 Branch points for kernel equation h{x,y) = 0 

For the censored random walk, we consider the kernel equation h{x, y) = 0 defined by the kernel function 
h{x, y) in the fundamental form (13.11) . Write ah{x, y) as a quadratic form in y with coefficients that are 
polynomials in x: 


ah{x, y) = (A 2 x)y^ + Aix^ - (Ai + A 2 + Ai + y 2 )x + yi 
For a fixed x, h{x, y) = 0 has two solutions 

-6(x) ± v^A(x) 


y + hx- 


(4.1) 


Y±{x) = 


2 A 2 X 


where h{x) = ab{x) = Aix^ — (Ai + A 2 + Ai + P' 2 )x + yi and A(x) = 6+(x)6_(x) with 
6+(x) = b{x) + 2x-\/A 2 A 2 = {x - l)(Aix - pLi) - {\[^2 - \^)'^x, 


6_(x) = b{x) - 2x\J A 2 A 2 = [x- l)(Aix - fii) - {^1^2 + y/^)'^x. 

Denote the branch points by Xi,i = 1, 2,3,4, which are the zeros of A(x), then we have 

6+(x) = Ai(x — X 2 )(x — X 3 ) and 6_(x) = Ai(x — xi)(x — X 4 ), 


(4.2) 


where 


0 < xi < X 2 < 1 < fj-i/Xi < xs < X 4 < +00 


and 


3:3 = 


(Ai + fii) + (\/^ — — \ (Ai + Ai) + (\/a2 — ~ 


2Ai 


(4.3) 


is a candidate of the dominant singularity of 7r^^^(x). 
Consider the following cut planes: 


Cx — Cx [^8)^4]) 

Cy = Cy [2/3 ) 2/4]) 

Cx = Cx- [X 3 ,X 4 ] U [xi,X 2 ], 

Cy = Cy - [2/3,224] U [2/1, 2/2], 

where and Cy are the complex planes of x and y, respectively. In the cut plane C^,, define the two 
branches of (F(x) by 


yo(a^) = F_(x) and yi(x) = F+(x) if |F_(x)| < |y+(x)|, 
^ 0 ( 2 :) = Y+{x) and Fi(x) = F_(x) if |y_(x)| > |y+(x)|. 


Symmetrically, when x and y are interchanged, we also have branch points yi,i = 1, 2,3,4, satisfying 


0 < 2/1 < 2/2 < 1 < 2/3 < 2/4 < +00 

as well as the two branches XQ{y) and Xi{y) defined in a similar fashion. 

Detailed properties of the branches Fo(2^) and Fi(x) {XQ{y) and Xi{y)) are needed in the asymptotic 
analysis for functions 7r|^^(x) and 7r|^^(x) {TT^\y) and TTY\y)), which are presented in the following two 
lemmas. 
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Lemma 4.1 The functions Yi(x), i = 0,1, are meromorphic in the cut plane Tlx- In addition, 

(i) Ib(a^) has one zero and no poles and Yi{x) has two poles and no zeros. Hence, Ib(a^) is analytic in 

Cx- 

(ii) l^o( x)! < |Li(x)| in the whole cut complex plane Cx- |5^o(2;)| = l^i(3j)l takes place only on the cuts. 

(iii) |lb( 3 ;)| <1 if\x\ = l,xjl^ 1, and Lb(l) = min ^ 1 ) < 1 - 

(iv) For all x G Cx, |Lb(x)| < and |Li(x)| > ■\J^- 

(v) If X & [xi,X 2 \, then |lo(x)| = and Xo(yo(®)) = x. 

V 

Moreover, 

(vi) 0 < lo(3;) < 1 for 1 < X < (recall that Ai < fii). 

Ai 

Parallel results for Xi{y), i = 0,1 can he stated as well. 

Proof. See [7], [I6] and m for proofs of (i)-(v). Here we only detail the proof to (vi). 

For 1 < X < let = 0 which leads to 

Aix + y = 0, 



where b{y) = \ 2 y‘^ - (Ai + A 2 + /li 


-r P2)y^ 1^2- 


on 


Ai Ai 

Ai’ Ai. 


, therefore, 2yJXifii < Aix + ^ < Ai + /ii. For y < 0, the inequalities 


/ b{y) + 2y\/Ai/li > 0 

\b{y) + (Ai + fii)y < 0 

have no solutions. For y > 0, solve the following inequalities 

/ b{y) +2y^Ai/li < 0 
\b{y) + (Ai + /ii)y > 0 

to have y 2 < y < ruin (l, or max (l, <y<yz- This means that for 1 < x < y 2 < Yq(x) < 1. 
^ A2 ' A2 Ai 

□ 

Lemma 4.2 We present more properties about Fo(x) and Xo(y) below: 

(i) ///i2 > A2, then 0 < To(3;) < 1 for x G (l, ^), and 1 < Yq{x) < for x G (^jXa). Specially, 

yo(f) = l andPoCxa) = > 1. 

(ii) If (12 < -^2) then 0 < Yq{x) < 1 for x G (1,X3). Also, Fo(A-) = A' < 1 and 10(2^3) = \fW < 1- 

(iii) If fi 2 = X 2 , then X 3 = 0 < Yo{x) < 1 for x £ ( 1 , and yo(l) = ^o(|^) = 1- 
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Similarly, 

(i’) If fi 2 > A 2 , then 0 < Xo{y) < I for y e (l, and 1 < Xo(y) < for y G (^,^ 3 )- Specially, 

Xo(g) = l anc?Xo(y3) = y/f > 1 . 

(ii’) If jl 2 < A 2 , i/ien 1 < Xo(y) < for y G ( 1 ,^ 3 )- Xo(l) = 1 and Xo{y 2 ,) = > 1. 


Proof. Based on equations (14.111 (14.2p . it is easy to know that h{x) < 0 for x G (1,X3), so the branch 
Yq{x) should take P-(x) = Solving the inequalities Po(a:^) > 1 and ^ 0 ( 2 :) < 1 with 

2A2X 

X G ( 1 ,X 3 ), we obtain, after some simple calculations, the results in (i)-(iii) of the lemma, (i’) and (ii’) 
can be proved in the same way. □ 


4.2 Poles of 7r^^^(x) 

Since the censored random walk is a standard walk in the quarter plane, literature results can now be 
applied to the analysis of the dominant singularity of 7 r^^^(x). Therefore, besides the branch point X 3 , 
given in ()4.3p . the other candidate for the dominant singularity is a pole of function 7 r|^^^(x). In the 
following, we rehne literature results, which lead to an explicit characterization of both the dominant 
pole and the regions for different exact tail asymptotic properties. 

Theorem 4.1 (Theorem 4.4 in |17] 1 If Xp is the pole of tt^\x) with the smallest modulus in {l,x^\, 
then Xp is a zero of hi{x,YQ{x)) orY^^Xp) is a zero of h 2 {XQ{y),y). In the latter case, |yo(3;p)| > 1- On 
the other hand, if Xp is the zero of hi{x,YQ{x)) orYQ{xp) with |lo(a^p)| > I is a zero of h 2 {XQ{y),y) with 
the smallest modulus in ( 1 ,X 3 ], then Xp is the pole of 7r^\x) in (IjXs]. Moreover, Xp is real. Parallel 
results can be easily stated for TT^\y). 

For the retrial queue model with two input streams and two orbits studied in this paper, we show 
in the following that the pole of 7 rj^^(x) (respectively 7r^\y)) can only be the zero of hi{x,YQ{x)) 
(respectively h 2 {Xo{y), y)). 

First, we discuss properties of the pole of 'k^^\x) in interval (IjXs). For convenience, let x* be 
the unique zero in (IjXs] of the function hi{x,YQ{x)) if such a zero exists, otherwise let x* = +00 (in 
this case, obviously x* can never be the dominant singularity since X 3 < + 00 ). Instead of directly 
considering the equation hi{x,YQ{x)) = 0 , we consider the product of two functions hi{x,YQ{x)) and 
hi{x,Yi{x)), which is a polynomial: 


hi{x,YQ{x))hi{x,Yi{x)) 


h2 

o:(A + yiY 


{x 




where 


g{x) = AAi(A + pLi)x^ + \yi{\ + yi- pi)x - yy.'l. 


Since A(Ai +^i) < /r/ii, it is easily to check that g{0) < 0 and g{l) < 0. Hence g{x) = 0 has one positive 
zero x+ and one negative zero x_. Especially, x+ > 1. The expressions of the two zeros are given as 


—A//i(A + gi — y) + y^[A/xi(A Y yi — /r)]^ + 4AAi(A + yi)yy\ 

2AAi(A + yi) 


(4.4) 


-\yi{\ + yi- y)- v^[A^i(A + yi- y)Y + 4AAi(A + yi)yyi 

2AAi(A + yf) 
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Since either pjj or p\ j is not X-shaped (refer to |17] for details) in this censored random walk, 
Theorem 4.5 in Li and Zhao m guarantees that the candidate zero of hi[x,YQ{x)) can only be x+. 
Solving /ii(x+,y) = 0, and then from (13.3p we get 


y = Y{x+) = 1 - 


[Ai(A + pi)x+ - piHi]{x+ - 1) 

A2(A + /ii)x+ 


(4.5) 


where T(x+) is either To( 2 :+) or Yi{x+). On the other hand, ^ ^ g( Ai(a+Vi) ) ~ > 0, 

hence, 1 < x+ < means Y{x+) > 1. Furthermore, to check whether or not x+ is the pole 

of 7 r|^^(x), we will carry out a discussion under the condition fi 2 > A 2 and /i 2 < A 2 , respectively, in the 
following lemma. 

Lemma 4.3 1. When (I 2 > X 2 , the value of x* depends on the value of x+: 


(a) For € (1, ^ , we have x* = + 00 ; 

\ J 


(b) For min (x 3 , Ai(a+Vi) ) j’ have x* = x+ < X 3 ifY{x+) < 

otherwise; 

(c) For x+ = X 3 < a^(a+Vi) ’ have Y{x+) = and x* = x+ = X 3 ; 

(d) For X 3 < x+ < we have x* = + 00 . 


and X* 
A 2 


= +00 


2. When P 2 < X 2 , we have x* = + 00 . 
Proof. For the case 1 X 2 > A 2 , if x+ G 


it leads to 0 < Yq{x) < 1 from Lemma For 


the case jl 2 < A 2 , if x+ G (1,X3], it leads to Yq{x+) < 1 /^ ^ 1 from Lemma [4.11 The both cases 

V ^2 

contradict T(x+) > 1. Hence, we can conclude that hi{x,YQ{x)) has no zero on [l,+oo) (x+ is the zero 
of hi{x,Yi{x))). Therefore, x* = + 00 . Other conclusions are easy to make. □ 


Remark 4.1 It is worthwhile to notice that: (i) If there does not exist a pole in ( 1 ,X 3 ], then X 3 is the 
dominant singularity of 7 rj^^(x). Therefore, for the purpose of dominant singularity, we do not need to 
consider case 1 (d) in Lemma 031 (h) The right-hand expression in (14.5p can be either yo(3:-i-) or Y’i(x+). 
(hi) It is possible that both x+ and x_ are zeros of hi{x,Yi{x)). In this case, hi{x,YQ{x)) = 0 has no 
solution. 

Next, we show /i 2 (Xo(g)) 2 /) has no zeros. Based on Theorem 14.11 if the pole in (IjXs] of 7 r[^^(x) is 
not X*, then it is denoted by xi. For convenience, define y* to be the unique zero of h 2 {XQ{y),y) in 
( 1 , 2 / 3 ] if such a zero exists, otherwise let y* = -|-oo. Following the same idea as above, we have 


h2{Xo{y),y)h2{Xi{y),y) = ^ —^(y - l)/(y), 

a(A p.2r 

where 

f{y) = AA2(A + fi2)y‘^ + Ay2(A + 112- ti)y - Fhl 

has two zeros: y_ < 0 and y+ > 1. Solving /i 2 (x,y+) = 0, and then from ()3.4I) we get 


X = X{y+) = 1 - 


[A2(A + y2)y+ - yF2]{y+ - 1) 

Xi{X + p2)y+ 
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where X{yj^) is either Xo(y+) or Xi{y^). 

Using a similar argument, parallel results to Lemma H131 -1 can be obtained. Since A 2 (A + 112 ) < fJ'IJ -2 

and /( ;v 2 (a+/. 2 ) ) = 1 < y+ < a 2 (a+/. 2 ) - ^(y+) > 1- Next, we claim 

that xi cannot exist. 

If h 2 {Xo{y),y) has a zero y* in ( 1 ,^ 3 ], then y* = y+ and 1 < X{y+) = Xo{y*) < From 


Theorem 4.7 in m we know xi 


Xi{y*). Hence, xi G 




and 0 < 10 ( 5 ^ 1 ) < 1 from Lemma 


14.21 Obviously, it contradicts to that To(5i) is a pole of 7r2^^(y). Therefore, xi cannot exist. 

Based on the above discussion, we are ready to summarize the detailed properties on the location 
of the dominant singularity. For convenience, we introduce the following three conditions: 


Condition 1. 
Condition 2. 


fi 2 > A 2 , x+ G (g, min (x 3 , A,(l+Vi) )) ^(^+) < 75 

fi 2 > X 2 and x+ = X 3 G (g, a,(a+Vi) )- 


Condition 3. One of the following three: (a) /i 2 < A 2 ; (b) jl 2 > X 2 and x+ G (1, g]; and (c) fl 2 > A 2 , 

1 +e (fe.min (X 3 , j^)), and r(l+) > 

Lemma 4.4 Case 1: Under Condition 1, the dominant singularity Xdom = x* = Xj^ < X 3 , which is a 
pole. 


Case 2: Under Condition 2, the dominant singularity x^om = X 3 = x* 
point and a pole. 


which is both a branch 


Case 3: Under Condition 3, the dominant singularity Xdom = x^ < x* = + 00 , which is a branch point. 

Remark 4.2 One should notice that the above lemma is a refinement of the literature result for a gen¬ 
eral random walk in the quarter plane. It provides explicit conditions (in terms of system parameters), 
under which the dominant singularity Xdom of vrg (x) (also explicitly expressed) is either Xdom = or 
Xdom = since all the branch point X 3 , the pole and U(x+) are explicitly expressed in (14.3p . ()4.4p 
and (14.5p . respectively. 


5 Asymptotic properties of 7T^\x) at its dominant singularity 

Once again, in this section, we only provide detailed analysis for the function 7 rg(x). Due to symmetry, 
parallel results for Tr^\y) can be easily stated and similarly proved. In the previous section, we proved 
that either X 3 or is the dominant singularity of 7 rg(x). In this section, we prove (in Theorem 15.ip 
that there exist three types of detailed asymptotic properties as x approaches to the dominant singularity 
Xdom of 7 rj^^(x), depending on Xdom = x+ < xs or Xdom = X 3 < x+ or Xdom = x+ = X 3 respectively. 

For simplicity in the following discussion, especially for the case of Xdom = X 3 , we write 

Yo{x) = p{x) + q{x)./l - — , (5.1) 

V ^dom 


hi{x,YQ{x)) =pi{x) +qi{x). I - 


^dom 
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doin') ^o{^) 


1 - 


X 


^dom 


p*ix)-q{x)Jl- 


^dom 


/il(x,yo(x)) ^(^dom)) 


1 - 


^dom 


Pi{x) + qi{x). I - 


^dom 


where 


b{x) 1 A{x) 


, if XdQui X 3 , 


, , h(x)ai{x) , , , , , / N / N 

Pi[x) = --hftiix), qi[x) = ai{x)q{x), 

2 A 2 X 


*/ \ ^Pixdom) p{x))X(iom 1 \ iPl{x) Pl{Xdom))Xdom 

p [X) = - and Pi{x) = -. 

^dom ^ ^dom ^ 

Theorem 5.1 The behaviour of tt^\x) at the dominant singularity is given as 
(i) Ifxdom = X* = x+ < X3, then 


lim (^1 - —= Ci^o, 
x^x+ y x+J 


where 


Ci,o = 


(A + /ii)v^A(x+) /i2(x+,To(x+))4^^(yo(a:+)) + /io(a:+,yo(a;+))7ro,o(l) 


/i 2 AAix+(x+ — l)(x+ — X-) 


(ii) If Xdom = X3 = X* = x+, then 


lim - x/xdom'^\^\x) = C 2 , 0 , 

^ ^^dom 


where 


C 2 ,Q — — X 


2 h2(^Xdomi^oiXdom))^2 (^(^rfom)) “1“ ^o(®dom) ^(^dom))^0,o(f) 


\/Xdom {xdom Xi ) {xdom ^2) (®4 Xdom) 


(iii) If Xdom = X3 < x* = + 00 , then 


lim [x) = Cs^, 




where (x) is the derivative of {x) and 


Csfi — — 


qix3) d 

2x3 dy 


h2{x3,y)-K)^\y) + /to(x3,?/)7ro,o(l) 
hi{x3,y) 


y=yo{^3) 
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Proof, (i) If Xdom = X* = x+ < X3, then Xdom is a simple pole of 7rj^^(x). Based on the analysis in 
m , we can rewrite 


7rf)(a;) = - 


h2ix,Yo{x))7r^2\^oix)) + ho(a:,yo(3:))7ro,o(l) hi(a;,Yi(x)) 


hi{x,YQ{x))hi{x,Yi{x)) 


h2{x,Yo{x))Tr^^\Yo{x)) +/io(x,lo(a^))7ro,o(l) hi(x,Yi(x)) 


J(aw(^ - l)g{x) 


h2{x,Yo{x))Tr^^\Yo{x)) +/io(x,lo(a^))7ro,o(l) hi(x,Yi(x)) 


- l)AAi(x - x-){x - x+) 


It follows that 


x— 


X+ 


lim I 1 -) vrj {x) = Ci 


(ii) If Xdom = X3 = X* = x+, then hi{xdom,YQ{xdom)) = 0 . In this case, we can rewrite 7 rj^^(x) as 

„(i)/ ^ _ -^2{x,Yq{x))'k^2\^o{x)) - /io(x,lo(a;))7ro,o(l) 

— 


V^l - xjxdom a/1 - xjxdamVX (^) + (®) 


It follows that 


lllTL v ^f ^dom'^X — 7 \ 7 \ 

(X\ (^dom)^(^cZom ) 


^ ^^dom 

(hi) If Xdom = X3 < X*, let 

T{x,y) = 


— C 2 , 0 - 


-h2{x,Yo{x))Tr^2'^ (Yo{x)) - /io(x, lo(a:))vro,o(l) 


hi{x,YQ{x)) 


Then the derivative of (x) is given by 


^(1), , dT , dTdYoix) 


with 


dYo{x) 


= p'ix) + q'{x)y/l - x/xdom “ 


q{x) 


‘^Xdom X!Xdom 

where p(x) and g(x) are defined in equation (j 5 . 1 |) . Since it is obvious that lim — xjx^~ 
“^ST’ \/l - xjxz^ = 0 and ^ is continuous at (xs, 10(2:3)), so, 


21 —>-X 3 


! - /Q\ 

lim ^J\-xIx3'k{’{x) = - 

2 X3 oy 


q{x 3 ) dT 


= C 3 . 


0 - 


(xiXoi^s)) 


□ 
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6 Tail asymptotic properties in stationary probabilities 


Exact tail asymptotic properties in stationary probabilities are obtained directly from the corresponding 
asymptotic properties of the unknown generating function by applying the following Tauberian-like 
theorem. This theorem is originated from Bender [2], and more complete versions can be found in 
Flajolet and Sedgewick [8], which include the following theorem as a special case. 


Theorem 6.1 (Tauberian-like theorem for single singularity) Let A{z) = analytic 

at zero with R the radius of convergence. Suppose that R is a singularity of A{z) that can he continued 
to a A-domain at R. If for a real number /3 ^ {0, —1, —2,...}, 

lim (1 — zlR)^A{z) = g, 


where g is a non-zero constant. Then, 




where r(/ 3 ) is the value of Gamma function at f3, and On ~ bn is equivalent to lim„a„,/ 6 „ = 1. 


The Tauberian-like theorem claims that the tail behaviour in the sequence of the coefficients in the 
Taylor expansion of the analytic function corresponds to the asymptotic property of the function at its 
dominant singularity. In the following subsections, we show how to apply Theorem 16.11 to characterize 
the tail behaviour in the joint probabilities T^m,n{k) for a hxed number n of customers in orbit 2 . 
Specifically, in subsection 1 , we provide a characterization for tail asympotics, when the server is busy, 
in the sequence of: ( 1 ) boundary probabilities 7rm,o(l); ( 2 ) marginal probabilities iTm = '^m,n(l); 

(3) joint probabilities 7 rm^„(l) for a fixed n > 0 (along the direction of queue one). While in subsection 2, 
when the server is idle, we provide a characterization for tail asympotics in 7rm,n(0) for a hxed n and 
for the marginal distribution TTm^ = 

Remark 6.1 By symmetry, tail behaviour in '7rm,n(l) and 7 rm^„( 0 ) for a hxed number m of customers 
in orbit 1 (and also in the marginal distributions for the second queue length when the server is busy 
and idle, respectively) can be easily stated and similarly proved. 


6.1 Exact tail asymptotics when the server is busy 

First, we consider the sequence 7rm,o(l) of the boundary probabilities. When the second queue is 
empty and the server is busy, the exact tail asymptotic behaviour of the stationary probability sequence 
7rm,o(l) along the increasing direction of the hrst queue is a direct consequence of the characterization 
of the asymptotic property for the function (x) in Theorem 15.11 and the Tauberian-like theorem 
(Theorem 16.Ih . 

Theorem 6.2 For a stable retrial queue with two input streams and two orbits studied in this paper, 
when m is large, we have three types of tail asymptotic properties for the boundary probabilities 7rm,o(l)- 

Type 1: (Exact geometric decay) Under Condition 1, 

( 1 

7rm,o(l) ~ Cifl ( — 1 I m > 1; 

Type 2: (Geometric decay with prefactor Under Condition 2, 

,,, C2,0 -l( 1 

7rm,o(l) ~ 2 - , m > 1; 

Y TT y ^dom J 
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Type 3: (Geometric decay with prefactor m 


Under Condition 3, 


7rm,o(l) ~ 



3 

2 




m>l. 


Here, constants Ci^ (i = l,2,3j are given in Theorem \5.1\ 

Remark 6.2 One may notice that in Type 3, the power of the decay rate is m — 2 instead of m — 1 
since the Tauberian-like theorem is applied to the derivative of the function. 

For characterizing the asymptotic behaviour of the marginal probability 7rm,n(l)) we 

compute 7 r(^)(x, 1 ), 


7r^^^(a:, 1) 


hi{x, l)7rf^(a:) + h 2 {x, l)7r^^^(l) + ho{x, l)7ro,o(l) 

—h{x, 1) 

[Ai(A + /ii)x - pi] 7ri^^(x) + Ai4^^(l) + Ai7ro,o(l) 

a [Ai(A + pi)a: - pi] 7r[^^(3;) + Ai(A + pi)(7r^^\l) + 7ro,o(l)) 
A + W 


If A 2 7 ^ P 2 ) it follows from (I4.2p that we have pi/Ai < 3 : 3 . Therefore, from Lemma 14.41 we can claim 
that 1 < pi/Ai < min(x*,X 3 ) is always true. Obviously, pi/Ai is the dominant singularity of 7 ri^i(x, 1), 
which is a simple pole. If A 2 = P 2 ) then from Lemma [4.2I liiii. we have X 3 = pi/Ai. Again, according to 
LemmaSUl the dominant singularity of 'k^^\x) is X 3 = pi/Ai < x* = + 00 . Notice that lim 2 ;_>.a ;3 7 rj;^^(x) 
is hnite. Therefore, the Tauberian-like theorem can be still applied. 


Theorem 6.3 (i) 

lim_ (l - ^ ^ 1 ) = Cm 

\ pi/Ai/ 

where 


Cm = +7ro,o(l)); (6.1) 

A + pi Pi ' 

and (a) The marginal probabilities has an exact geometric decay with decay rate equal to x^om = 

pi/Ai.- 


vr 


( 1 ) 


~ Cm 



Remark 6.3 It should be noticed that one may consider +'^m,n{0)) the usual marginal 

distribution of the first queue. Its tail asymptotic property can be easily obtained since the property 
for and 7rm,o(l) have been studied, and the property for = X^^i^m,n( 0 ) and 7 rm,o( 0 ) can be 
similarly obtained. 

Next, the exact tail asymptotic behaviour for joint probabilities can be obtained from the recursive 
relationship of the generating functions (pn{x), defined by 


00 

lfn{x)='^7rm,n{i)x""~^, U > 0. 
m=l 
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It is clear that i^o{x) = '7r[^^(a;). From the balance equations of the censored random walk, we can obtain 


c{x)(fi{x) + 6i(x)<y9o(x) =ao(x), 
c{x)ip 2 {x) + h{x)ipi{x) + ai{x)ipQ{x) =al{x), 
c{x)(pn+i{x) + b{x)ipn{x) + a{x)ipn-iix) = al{x), n>2, 


( 6 . 2 ) 

(6.3) 

(6.4) 


where 


Rewrite (16.4p as 


aQ{x) = -C2(x)7ro,i - 6o(3:)7ro,o, 

«i(a;) = -C2(a;)7ro,2 - b2{x)7ro^i - ao(x)7ro,o, 

o* (x) = -C2(a;)7ro,n+i - &2(a;)7ro,n - a2(x)7ro,n-i, n > 2. 


¥^n+l (x) 


-b(x)(pn(x) - a(x)(fn-i(x) + a*(x) 
c(x) 


n > 2, 


and note that c(x) = po,-ix. Hence, we established the fact that ipn{x) has the same singularities as 
<y9o(x) since that the zero of c{x) is not a pole of ipn{x) for all n > 0. 

By adopting Theorem 7.1 and Lemma 7.2 in m directly, we define 


^iixdom) — 


i^dom) ^ 
7 

C\X(lom ) 


i = 1,2,3 


and 


Bsixs) = 


Pijxs) 

c(x 3 ) 


^3,0, 


then we can conclude the results in the following theorem. 

Theorem 6.4 Corresponding to the three types in Theorem \5 . R when m is large, we have the following 
tail asymptotic properties for the joint probabilities 7rm,n(l) for a fixed n: 

Type 1: (Exact geometric decay) 


'^m^n (1) ~ Ai{x+) 


n—1 


yiix+) 

Type 2: (Geometric decay with prefactor 


1 

x+ 


m—1 


, n>l; 


-^ 2 (^dom) ( 1 


n—1 


w(i)~ ^ 


Type 3: (Geometric decay with prefactor m 


m 2 


^dom 


m—1 


, n > 1; 


^m,n(l) 


[A3 {x3 ) + {n - 1)B3 {x3 )] f 1 


vr 


\yi{x3) 


n—1 


m 2 I — 

.^3 


771—2 


n > 1. 


6.2 Exact tail asymptotics when the server is idle 

Having known the exact tail asymptotic properties of the boundary, marginal and joint distributions 
for I(t) = 1 (or the server is busy), we can now study the tail asymptotic properties for I{t) = 0 (or 
the server is idle) based on the relationship given in (|2.2p . 

Setting y = 0 in (12.2p leads to 

(A + /ri)pW(a:,0) =/rP(i)(x,0) +//iPW(0,0), (6.5) 
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which means that and have the same asymptotic property. 

Similarly, setting y = 1 in (|2.2I) leads to 


1 ) = 1 ) + H2P^^Hx,0) + ( 0 , 1 ). ( 6 . 6 ) 

Substituting (j6.5p into ()6.6p gives 

aP^°\x,l) = p(^^(x,0) + +yiP(°)(0,l). 

A + y-i A + Hi 

Since the asymptotic property at the dominant singularity of P^®^ {x,l) is dominated by the asymptotic 
property of the function yP*^^)(x, 1), P^^\x, 1) and P^^^(a;, 1) have the same asymptotic property. Based 
on the above, we have the following conclusion: 

Theorem 6.5 Assume that the retrial queue with two input streams and two orbits is stable. 

(i) For large m, we have three types of tail asymptotic properties for the boundary probabilities 7rm,o(0) 
correspondingly. 

Type 1: (Exact geometric decay) 


7rm,o(0) ~ , ^ Ci ,0 ( , m > 1; 


X + Hi ’ \3:^+ 

Type 2: (Geometric decay with prefactor m~^P) 


m—1 


'^m,o( 0 ) 


H C2,o -if 1 

rm 2 ' 


A -|- Hi ypK \^dom 

Type 3: (Geometric decay with prefactor m~^P) 


m—1 


m > 1; 


7rm,o(0) 


H C's.o -I ( ^ 

- - -^m 2 — 

A + /U1 Vvr Vx3 


m—2 


m > 1. 


Here, constants Ci^ {i = 1,2,3) are given in Theorem \5.1[ 

(ii) The tail asymptotic property of the marginal distribution '^rn,ni0) is determined by 


m—1 


a ypi J 

where Cm is provided by dSU). 

We finally characterize the tail asymptotic behaviour for the joint probabilities TTm,n{^) for a fixed 
n > 0. Define the generating function 

CX) 

Gn\x) = ^ 7rm,n{k)x^, k = 0,1 n > 1. 

m=0 

Referring to equation (14) in [3], we have 

aG''^\x) - HGff\x) = //i7ro,„(0), 
which obviously leads to the following theorem. 
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Theorem 6.6 Corresponding to the three types in Theorem \5 . R when m is large, we have the following 
tail asymptotic properties of the joint probabilities for a fixed n: 

Type 1: (Exact geometric decay) Under Condition 1, 


R 


vrm,n(0) ~ -Ai{x+) 


a 


1 


Yi{x+) 


n—1 


1 

x+ 


m—1 


, re > 1; 


Type 2: (Geometric decay with prefactor rre Under Condition 2, 

1 


R A2(xdoim) / 1 _ 

vrm,n(0)~-W- n' 

OL -y/TT 


^dom 


n > 1; 


Type 3: (Geometric decay with prefactor rre Under Condition 3, 


'^m,n ( 0 ) 


uiAsixs) + {n-l)B3{x3)] f 1 


a 


IT 


\Yi{x3) 


n—1 


m—2 


rre 2 [ — ] ^ re > 1. 


7 Concluding remarks 

In this paper, we considered the exact tail asymptotic behaviours of a retrial queue with two input 
streams and two orbits. Partitioned according to the two states of the server, this model is formulated 
as a random walk in the quarter plane whose transition probabilities are modulated by a two-state 
Markov chain (idle or busy). Our work is a revisit of the same model studied in [3]. While in [3], 
the study is based on the solution to a BVP, we employed a different method, the kernel method. 
The advantage of using this method mainly relies on the fact that there is no need to have a full 
determination of the unknown generating function. Instead, we only need the location of the dominant 
singularity of the unknown function and the asymptotic property at its dominant singularity. By this 
method, tail asymptotic properties in stationary probabilities for the model are obtained when the first 
queue size is large. Due to symmetry, it is not difficult to state and (similarly) prove parallel exact tail 
asymptotic properties when the second queue size is large. In addition, exact tail asymptotic results for 
other probability sequences formed from the joint stationary probabilities can also be considered. For 
example, we can consider the total number of customers in the system as follows: let 

— 'y ^ 'Xm,n 

m,n\ 

m-\-n=T 


and we compute 'k^^\x,x), according to m-. 

(Ax - + {Xx - yf^)4^^(x) + A7ro,o(l) 


( 1 ) 


{x,x) = -- 


x{Xx — 


A 1 +A 2 ' 


Then, the dominant singularity is determined by comparing x = {fii+fi 2 )/X to the dominant singularities 
of 7r(^^(3:) and 7r2^^(x), and therefore the asymptotic property at its dominant singularity is determined. 
The exact tail asymptotic property is a consequence of the Tauberian-like theorem. 

This paper used a censored chain to convert the matrix-form fundamental form into a usual (scalar) 
fundamental form. It is not always feasible to do this conversion since explicit expressions might not 
exist for the censored chain. A general method is to solve the matrix-form fundamental form to have a 
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relationship between generating functions for different states of the modulated chain. For example, the 
censored chain to the idle state does not have an explicit expression for its transition matrix. However, 
in terms of the relationship in ( 12 . 2 |) and (j2.3|) obtained by solving the matrix-form fundamental form, 
we can have the following functional equation: 

R{x,y)P^°'>{x,y) = A{x,y)P^'^\x,0) + B{x,y)P^°\o,y), |x| < l,\y\ < 1, 

with 


R{x, y) = Ai(l - x)xy A 2 (l - y)xy - yi{l - x)y - /l 2 (l - y)x, 
Mx, y) = [(1 - y){My - /i) + Ai(i - x)y] y 2 x, 

B{x, y) = [(1 - x){\ix - y) + A 2 (l - y)x] yiy, 


which is equivalent to: 

^ Mx,y)-Rix,y) ( 0 ), ^ ^ B{x,y) - R{x,y) (q) A{x,y) + B{x,y) - R{x,y) 

R{x,y)P >{x,y) = -vr) '(xjH-vr^ '(y)H-vro,o(0). 

y X xy 

After some calculations, the above equation also can be written as 

-U^\x,y)-n^^\x,y) = hf\x,y)T:f\x) + 4°^ (x, -h Ji§\x,y)-Kofl{^), 

where 


(x, y) = [Aix A 22 / ilix~^ -h il 2 y~^ -{\ + (jLi+ il 2 )]xy, 
hf\x,y) = [Ai(A -h yLi)x -h A 2 (A -h yLi)y + fiix~^ - A(A yi) - jli]x, 
hf\x,y) = [Ai(A -h y 2 )x -h A 2 (A -h y 2 )y + y 2 y~^ - A(A -h ^ 2 ) - y 2 ]y, 

(x, y) = AAix -I- X\ 2 y - A^. 

The above functional equation is the fundamental form corresponding a random walk defined by 
Pi,o = Ai, Po,i = A 2 , p-ifl = yi, po ,-i =/I 2 , Po,o = 1 — (A + Ai + A 2 ), 


Pifi = Ai(A-hpi), 

pS = A2(A-hpi), 

qi) 

P- 1,0 = 

Po,0 — ^ ~ + Pi) + PlJj 

pffl = Ai(A-hp2), 

Poj = X2(X + y2), 

~(2) 

P0,ll = P2, 

Po^o ~ ^ ~ + P 2 ) + /I 2 ], 


AO) _ AO) 

7*1,0 ~ 7*0,1 

— AA 2 , Po ,0 ~ 

1-A2. 


We now can apply the kernel method to the resulting fundamental form to obtain exact tail asymptotic 
properties for probabilities with an idle server. 

Finally, we emphasize that this work serves as an illustration of how the kernel method can be applied 
to random walks modulated by a finite-state Markov chain, which has a similar structure property to 
that the retrial queue model possesses. 
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